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SECTION  1  -  GENERAL  THEORY 


Because  of  the  intercontinental  ballistic  missile  and  satellite  programs, 
the  hypersonic  flow  around  a  body  of  revolution  has  been  investigated.  Recently 
the  wake  of  such  a  body  in  a  hypersonic  flow  becomes  one  of  the  most  interesting 
problems  among  aerodynamicists  due  to  the  desire  to  monitor  the  re-entry 
bodies.  The  thermal  radiation  emitted  by  the  hot  gas  left  in  the  trail  and  the 
electron  density  in  the  trail  are  the  observables  in  the  trails  of  such  bodies. 

In  order  to  determine  the  observables,  we  have  to  know,  in  detail,  the  flow 
field  of  the  wakes  behind  the  body,  which  depends  mainly  on  the  flow  in  the 
boundary  layer  of  the  body.  The  flow  field  in  the  boundary  layer  of  the  body  as 
well  as  that  in  the  wake  may  be  either  laminar  or  turbulent.  The  properties  of 
the  observables  depend  greatly  on  whether  the  flow  in  the  wake  is  laminar  or 
turbulent.  Hence  it  is  extremely  interesting  to  know  the  transition  point  or  the 
transition  region  of  the  flow  field  on  and  behind  a  body  of  revolution.  Even 
though  there  is  no  theoretical  analysis  available  to  predict  accurately  the  transi¬ 
tion  point,  the  stability  of  a  laminar  flow  has  a  close  relation  with  the  transition. 
In  this  report,  we  will  study  the  stability  of  axisymmetrical  flow  in  general, 
which  includes  both  the  stability  of  axisymmetrical  wake  or  jet  and  of  boundary 
layer  flow  on  a  slender  body  of  revolution. 

A  detailed  analysis  of  the  stability  of  a  hypersonic  flow  over  a  body  is 
very  complicated  because  not  only  is  the  flow  compressible,  but  the  composition 
of  the  air  may  vary  considerably  from  place  to  place.  An  accurate  analysis 
should  include  both  the  effect  of  compressibility  and  the  effect  of  th'  composi¬ 
tion  of  the  gas  in  the  flow  field.  However,  since  the  stability  of  the  axisym  ¬ 
metrical  flow  of  an  incompressible  fluid  has  not  yet  been  well  develope^^,  it 
would  be  useful  to  investigate  this  case  first  to  bring  out  the  essential  points  of 
stability  of  an  axisymmetrical  flow  before  an  extensive  program  of  the  stability 
of  an  axisymmetrical  flow  of  a  compressible  gas  with  the  influence  of  variable 
composition  is  carried  out. 
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The  problem  is  formulated  in  a  general  manner,  so  that  the  results  may 
be  applied  to  both  flow  problems  with  solid  boundaries,  such  as  boundary  layer 
flows  and  pipe  flows,  and  those  without  solid  boundaries,  such  as  wake  flows 
and  jet  flows. 

Since  we  consider  only  the  axisymmetrical  basic  flow,  the  general  dis¬ 
turbances  will  be  rotationally  symmetric  disturbances.  We  shall  formulate  the 
problem  of  the  stability  of  an  axisymmetric  flow  subjected  to  a  rotationally 
symmetric  disturbance  in  Section  2.  But  we  shall  show  in  Section  3  that  the 
axisymmetrical  flow  is  always  stable  with  respect  to  disturbance  of  the  tangential 
velocity  component.  As  a  result  we  need  only  to  study  the  stability  with  respect 
to  axisymmetrical  disturbances  in  detail  in  Sections  4  to  6. 

The  stability  equation  of  axisymmetrical  disturbances  is  investigated  and 
its  solution  for  large  Reynolds  number  is  found  for  arbitrary  basic  velocity 
profile  in  Section  4.  The  boundary  value  problems  of  the  stability  for  various 
axisymmetrical  flows  are  analysed  in  Section  5.  However,  the  detailed  numerical 
calculations  will  not  be  given  here  but  in  later  parts  of  this  report.  The  com¬ 
parison  of  the  stability  of  axisymmetrical  flow  with  that  of  two-dimensional  flow 
will  be  discussed  according  to  an  equivalent  principle  in  Section  6. 
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SECTION  2  -  GENERAL  FORMULATION  OF  THE  PROBLEM 


We  use  cylindrical  coordinates  r,  9,  z  with  velocity  components  u,  v,  w 
in  the  direction  of  r,  9,  z, respectively.  For  an  incompressible  and  viscous 
fluid,  the  equation  of  continuity  and  equations  of  motion  in  the  cylindrical  co¬ 
ordinates  are,  respectively, 
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For  an  axisym metrical  flow  subjected  to  rotationally  symmetric  disturb¬ 
ances,  all  the  variables  are  independent  of  the  angular  displacement.  Hence, 
^equations  (2.1)  to  (2.4)  reduce  to 
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p  is  the  density,  p  is  the  pressure  and  v  is  the  coefficient  of  kinematic  viscosity 
of  the  fluid.  Equations  (2.  5)  to  (2. 8)  are  the  basic  equations  for  our  stability 
analysis. 

We  consider  only  the  stability  of  parallel  or  nearly  parallel  flows  which 
include  flows  of  the  types  of  boundary  layer,  wake,  and  jet.  The  basic  velocity 
profile  of  these  parallel  and  nearly  parallel  flows  may  be  considered  as  follows: 


U  =  V  =  0  ,  W  =  W  (r)  (2.9) 

From  now  on,  all  quantities  are  to  be  non-dimensional  such  that  all  velocities 
are  expressed  in  terms  of  a  characteristic  velocity  (z),  all  lengths  are  ex¬ 
pressed  in  terms  of  a  characteristic  length  L  (z),  the  pressure  is  expressed  in 
2 

terms  of  d  V  ,  the  time  is  expressed  in  terms  of  L/U  and  the  Reynolds  number 
o  o 

of  our  problem  is 

L  U 

R  =  — ^  (2.10) 
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For  nearly  parallel  flow.  and  L  are  usually  functions  of  the  axial  distance  z. 
We  shall  discuss  the  choice  of  proper  values  for  and  L  when  we  investigate 
specific  flow  problems  in  the  later  parts  of  this  report. 

For  a  rotationally  symmetrical  disturbance,  all  the  variables  are  inde¬ 
pendent  of  the  angular  displacement  9  .  Hence,  the  velocity  components  and  the 
pressure  may  be  written  in  the  following  form; 


u  (r,  z,  t)  =  u  (r)  e.xp  [iQ  (z  -  ct)J 

V  (i^,  z,  t)  =  v  (r)  exp  [iQ  (z  -  ct)  ] 

w  (r,  z,  t)  =  W  (r)  +  w  (r)  exp  [ia  (z  -  ct)  ] 

p  (r,  z,  t)  =  P  (z)  +  p  (r)  exp  [ia  (z-ct)] 


(2.11) 


Where  W  (r)  and  P  (z)  are  respectively  the  velocity  and  pressure  of  the  basic 
flow;  u  =  u  (r),  v  =  v  (r),  w  =  w  (r)  and  p  =  p  (r)  are  respectively  the  perturbed 
velocity  components  and  pressure,  which  are  functions  of  r  only;  Q  is  the  wave 
number  which  is  real  and  positive  and  c  is  the  complex  wave  velocity,  i.e. , 

c  =  c  +  ic.  (2.12) 

r  1 

with  i  ==  vrr 


Substituting  equation  (2.11)  into  equations  (2.5)  to  (2.8)  and  neglecting  the 
higher  order  terms  of  the  perturbed  quantities,  we  obtain  the  following  equations 
for  the  perturbed  quantities 

(ru)’  +  i  o  r  w  =  0  (2. 13) 
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I  a  (W  -  c)  u  =  -  p*  +  ;;  (u**  +  —  -  a  u  -  (2.14) 

R  r  ^ 

r 
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iQ  (W-c)w  +  uW  =  -  iap+  -  (w"  +  —  -  a  w)  (2. 15) 

ia  (W-c)r  =  ^  (v"  +  ^  -a^r-^)  (2.16) 

r 

where  prime  refers  to  the  differentiation  with  respect  to  r. 

It  is  interesting  to  notice  that  equations  (2. 13)  to  (2. 16)  may  be  divided 
into  two  independent  groups : 

(i)  Equation  (2. 16)  is  independent  of  the  other  three  equations.  Since 
equation  (2. 16)  consists  of  the  unknown  v  (r)  only,  this  equation 
may  be  used  to  study  the  stability  of  an  axisymmetrical  flow  to 
the  tangential  velocity  disturbances. 

(ii)  Equations  (2. 13)  to  (2.  15)  are  coupled.  These  three  equations 
may  be  used  to  study  the  stability  of  an  axisymmetrical  flow  due 
to  some  axisymmetrical  disturbances  u,  w  and  p. 

We  may  reduce  equations  (2. 13)  to  (2. 15)  into  a  single  differential  equation 
for  the  function 


f  =  r  u 


(2.17) 


It  is  easy  to  show  that  the  function  f  is  proportional  to  the  perturbed  stream 
function  of  the  axisymmetrical  flow. 

Eliminating  the  pressure  p  (r)  from  equations  (2. 14)  and  (2. 15)  we  have 
2 
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We  may  express  u,  w  and  their  derivatives  with  respect  to  r  in  terms  of  f  and 
its  derivatives  in  r  by  the  help  of  equations  (2. 13)  and  (2, 17).  As  a  result, 
equation  (2. 18)  gives 


O  f»  Wf  i 

(W  -c)(f  -a  £  --)  -f  (W  *^)  =  -  ^ 

riv  2  42  2  3 

(f  -2a  £»»  +  a  f)  --  -  a  £•)  +  (f*»  --)  (2.19) 

L  J 

Equation  (2, 19)  is  the  fundamental  equation  for  the  stability  of  an  axisymmetrical 

flow  subjected  to  axisymmetrical  disturbances.  It  is  interesting  to  notice  that 

if  those  terms  containing  explicitly  the  factor  r  vanish,  equation  (2. 19)  reduces 

to  the  well  known  Or r -Somme rf eld  equation  for  the  stability  of  two-dimensional 
2 

parallel  flows. 

For  the  flow  field  between  r  =  r  and  r  =  r  ,  equation  (2. 19)  is  then  to  be 
solved  under  the  boundary  conditions 

f  (r^)  =  0  ,  r  (r^)  =  0  ,  f  (r^)  =  0  ,  F  (r^)  =  0  (2.20) 


Since  equation  (2. 19)  is  a  fourth  order  ordinary  differential  equation,  there 
exists  a  fundamental  system  of  four  solutions  of  equation  (2. 19)  which  are  analytic 
functions  of  the  variable  r  and  of  the  parameter  c,  oc  and  aR.  Let  the  four 
solutions  be  f  ,  f  ,  f  ,  and  f  .  The  conditions  (2.20)  gives  the  determinantal 
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Since  the  function  F(C  ,  a,  a  R)  is  an  entire  function  of  the  variables  c,  cl.  and 
aR,  we  may  solve  for  c  and  obtain 

c  =  c  (a,  R)  (2.22) 

Since  o  and  R  are  taken  to  be  real  and  positive,  we  have  four  equations  (2.21) 
and  (2.22)  =  c^  {a,  R);  c  .  =  c  .  (a,  R)  (2.23) 

The  curve  ci  =  o  gives  the  limit  of  stability. 

The  general  solutions  of  equation  (2. 19)  for  arbitrary  basic  velocity  pro¬ 
file  W(r)  are  difficult  to  obtain.  In  Section  4,  we  shall  find  the  asymptotic  solu¬ 
tions  for  the  limiting  cases  of  very  large  a  R. 
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SECTION  3  -  STABILITY  OF  AN  AXISYMMETRICAL  FLOW  WITH  RESPECT  TO 
THE  DISTURBANCES  OF  THE  TANGENTIAL  VELOCITY  COMPONENT  V 


The  function  v  of  the  perturbed  tangential  velocity  component  of  equation 
(2. 16)  should  be  considered  as  a  complex  function  of  r.  If  we  multiply  equation 
(2. 16)  by  rv  where  bar  denotes  the  complex  conjugate,  and  integrate  over  the 
interval  (r  ,  r  )  witli  respect  to  r,  we  obtain 
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(  v’  V  +  V*  r)  dr 
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where  the  boundary  conditions  v’  (r  )  =  v’  (r  )  =  0  are  used. 

1  ^ 

Substituting  equations  (2, 12)  and  (3.  2)  into  equation  (3. 1)  we  have 


(3.3) 


From  the  imaginary  part  of  equation  (3.3)  we  have 
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rdr 
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Equation  (3.4)  gives  the  velocity  of  propagation  of  this  tangential  velocity 
disturbance . 


From  the  real  part  of  equation  (3.3),  we  have 


1 

/ 


o  V  V  .  , 

2  —  —  )  dr 


(v*  v’  r  +  0!  vvr  +  r 


c. 

1 


1 

a  R 


vv  r  d  r 


(3.5) 


Since  the  integrals  in  equation  (3.5)  are  all  positive  quantities,  is  then  always 

negative.  As  a  result,  the  tangential  velocity  disturbances  will  always  be  damped 

out  for  all  axisymmetrical  flows  subjected  to  rotationally  symmetric  disturbances. 
3 

was  Lew  who  first  showed  that  this  result  holds  true  for  axisymmetrical  jet 
flow.  Since  the  flow  is  stable  with  respect  to  v,  we  need  only  to  consider  the 
axisymmetric  disturbances  governed  by  equation  (2. 19) 
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SECTION  4  -  APPROXIMATE  SOLUTIONS  OF  THE  STABILITY  EQUATION  OF 
THE  AXISYMMETRIC  FLOW  (2. 19) 


We  introduce  a  new  variable  ^  such  that 


5  =  1  -  V  '"■*> 

r 

c 


where  r  =  r  when  W  =  c. 
c 

Equation  (2. 19)  in  terms  of  ^  becomes 


(W  -  c)  Lf  -  (1  -  i)  W"f  =  -—hr  LLf 


where  prime  refers  to  derivative  with  respect  to  |  and 


(4.2) 


f  (r)  =  f  (4) 


L  =  a  -  4) 
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(4.3) 


It  is  interesting  to  notice  the  similarity  of  equation  (4. 2)  with  the  Orr- 
Sommerfeld  equation  of  two-dimensional  parallel  flow.  If  the  factor  (1  -  4)  of 
equation  (4.  2)  is  replaced  by  unity,  equation  (4. 2)  reduces  exactly  to  the  Orr- 
Sommerfeld  equation  of  the  stability  of  two-dimensional  flow.  Because  of  the 
similarity,  w«  may  obtain  the  solutions  of  equation  (4. 2)  by  the  same  method  used 
in  the  analysis  of  Orr-Sommerfeld  equations.  The  mathematical  properties  of 

f 

the  solutions  of  equation  (4.  2)  are  similar  to  those  of  Orr-Sommerfeld  equations.  ^ 
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'rhere  exists  a  fundamental  system  of  four  solutions  of  equation  (4.2)  which  are 
analytic  functions  of  the  variable  i  and  of  the  parameter  c,  and  being, 

in  fact,  entire  functions  of  these  parameters. 

We  may  find  the  solutions  of  equation  (4. 2)  by  convergent  series  in  terms 
of  a  small  parameter  e.  First  we  make  a  change  of  variable 

4  =  en  ,  f  (C)  =  0  (n)  (4.4) 

Equation  (4. 2)  becomes 
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2 

W»»  =  W”  +  W'"  (en)  +  (e-n)  + 

O  O  2! 


(4.6) 


subscript  o  refers  to  value  at  ^  =  0.  The  solution  is  then  obtained  in  the  form 


M€)  =  0(Ti)  =  +  6  0^^^  (n)  +  e^0^^^(n)+.  .  . 


(4.7) 
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and  the  differential  equations  for  the  approximations  of  successive  orders  can 
be  obtained  by  substituting  equations  (4.6)  and  (4.7)  into  (4.5)  and  equating  all 
the  coefficients  of  the  various  powers  of  e  to  zero.  From  equation  (4.5)  we  find 
that  the  proper  choice  of  the  parameter  c  is 

1 

■  3 

e  =  (a  R^)  (4.8) 


The  differential  equations  for  the  functions  (ti),  0^^^  (r|),etc. ,  are  as 


follows: 
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s'"  =  W  11  0^"^"''  =  L,  (0),  (ns  1) 
o  (n-1) 


(4.9) 


where  L,  _  (0)  is  a  linear  combination  of  0^^^  (ti),  0^^^  (ti),  ...  0^”  (r|)  and 

(n-1) 

their  derivatives.  In  particular 


W 

L  (0)  =  (W  --7^)Ti^0  +W'’  0  +  2i  (110^®^^'^+ 0^®^’") 

o  ^  o  Is  o  oo 


(4.10) 


(6). 


It  is  interesting  to  notice  that  the  zeroth  order  solutions  0  r|  are  exactly 

2 

the  same  as  those  of  the  two-dimensional  case,  i.e. , 


=  TI  ;  03^®)  =  dr]  /  dn  •  (z) 


-1-00  4-00 


0  =  1  ;  0  =  y*  dTi  f  dr\'  r?  (z) 
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(4.11) 
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where  z  =  j  (W  ri  and  are,  respectively,  the  first 

and  the  second  kind  of  Hankel  function  of  the  order  1/3 . 

The  solutions  of  the  higher  order  terms  ns  1  of  equation  (4.9) 

are  different  from  those  of  the  two-dimensional  case  because  the  functions 
L^n_i)  different  from  the  corresponding  functions  in  the  two-dimensional 

case.  The  higher  order  terms  are  given  by 


The  convergent  series  (4.7)  is  convergent  provided  e  is  restricted  so  that 
the  series  (4.6)  are  convergent. 

From  the  study  of  the  stability  of  two-dimensional  flow,  we  know  that  it 
is  usually  more  convenient  to  use  asymptotic  series  for  numerical  purposes, 
particularly  in  dealing  with  boundary  value  problems.  It  was  found  that  the  most 
convenient  way  is  to  replace  the  first  two  solutions  and  by  the  cor¬ 
responding  asymptotic  solutions  f^^^^  and  We  may  approximate  the  four 

fundamental  solutions  |f  ,  f  ,  f  ,  f  I  by  the  functions  if  f  0 
(0) )  1  1  2  3  4)  I  1  2  3 

0^  ^  in  the  boundary  value  problems.  The  asymptotic  solutions  of  f^  and  f^ 

may  be  obtained  by  developing  f  ( ^ )  in  powers  of  (“  We  put 


Substituting  equation  (4. 13)  into  equation  (4.2)  and  collecting  terms  of  the 
-1 

same  powers  of  >  we  have 

(W-c)  [(1-4)  f^®^"  f^^^]  -  (l-4)W"f^®^  =  0  (4.14) 

(W-c)  [(1-4)  f^"^  -  (1-4)  = 

^(n-l)"']  ^  ^ 

The  equations  of  the  asymptotic  solutions  (4. 14)  and  (4. 15)  are  similar  to  those 
of  the  two ‘-dimensional  case,  but  are  not  exactly  the  same  even  for  the  zeroth 
order  term  because  f^^^  corresponds  to  the  first  two  terms  of  the  convergent 
series  (4.7),  i.e.,  f^^^=  +  Since  0^^^  is  different  from  the  cor  > 

responding  function  in  the  two-dimensional  case,  f^^^  is  different  from  the  cor¬ 
responding  function  of  the  two-dimensional  case  too. 

There  are  two  solutions  of  the  type  of  (4. 13),  which  are  known  as  inviscid 
solutions,  and  which  can  be  obtained  by  developing  f^^^  in  powers  of  The 

two  particular  integrals  are 

f^<^)  =  (W-c)  [h^  (4)  +  ai^h2  (0  +  a^^h^  (C)  + •  .  •] 

(4. 16) 

f2^®>=(W-c)  [k^(0  +  ai^k3(4)  +  o;^^kg(4)+.  .  .] 

(4. 17) 
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where 


h  (|)  =  1 
o 


2n  +  2 


-  / 


(  (W  -  o)- 
a 


(W-c) 

2n^ _ L 


(1-4) 


d4 


d4 


(4. 18) 


nS  0 


S<i)=  /  ^ 

*  {  (W  -  cf 

cL 


*^211  +  3 


/ 


4  (W-c) 

a 


/  “an. 

f  — tt-t: - d4 


L^a 


(1-4) 


d4 


(4. 19) 


0 


where  4^  i^ay  be  any  fixed  point  but  4^  1- 

Near  the  axis  of  symmetry  4=1.  we  may  approximate  the  velocity  profile 
in  the  following  form : 


W  (r)  = 


B  r 


2n 


(4. 20) 


where  B  and  n  are  constants  to  be  determined  from  the  actual  velocity  pro¬ 
file.  Hence, 


W(4)  -  c  =  Br^^"  [l  -(1  -4)”] 
W”  =  -Br^^“n  [(n-l)  (1  -  4)""^j 


(4.21) 
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Equation  (4. 14)  for  the  velocity  profile  (4. 20)  becomes 


+  n  (n-1)  (1-4)"'^  ji)=  0 
(4. 22) 

where  tj)  =  f  for  the  velocity  profile  (4.  20).  Equation  (4. 22)  is  the  Pretsch's 
equation^  which  has  two  solutions 

=  il>,  (V  =  e  C  P,  (C) 

1  1  o  1 

f2^®^  =  02^^)  =  [P2  (^)  + +Ki/)i  (4.23) 

for  4  ^  4  ^  f 

a 

where  P  (4).  P„  (4)  A  for  n  =  1  to  5  are  given  in  Reference  1.  P^  ( 4) 

P^  ( 4 )  a*"®  power  series  of  4  •  The  constants  e^  and  K  are  determined  by  join¬ 
ing  the  solutions  f/^^  and  at  the  fixed  point  4=4  »  i.e. , 

11  a 


1  -(1-4) 


■  ,  ^ 
(i-4)-i 


~OL. 


4  P 

1 


<«a> 


=  f 


(0) 


=  W  (4  )  -C 
a 


o  o 

f2^®N4a)  =  0  (4-24) 

We  shall  use  the  functioas  f/^^  of  equations  (4. 16)  and  (4. 17)  for  4  s:  ^  >  _oo 

1  a 

and  the  functions  of  equatioa  (4. 23)  for  1  ^  4  ^  4  I)'®  inviscid  solutions 

1  a 

f  and  f  ,  respectively 
X  ^ 
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Having  found  the  two  particular  integrals  of  f^^\  we  can  obtain  the  higher 
approximations  by  quadratures.  However,  in  the  calculation  of  stability  of  an 
axisymmetrical  flow,  it  is  sufficient  to  use 


We  may  obtain  the  viscous  solutions  of  the  asymptotic  series  in  the  same  manner 
as  in  the  two-dimensional  case.  Since  we  are  not  going  to  use  these  asymptotic 
viscous  solutions,  we  shall  not  give  them  here. 
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SECTION  5  -  THE  BOUNDARY  VALUE  PROBLEM 


Having  obtained  the  fundamental  solutions,  we  may  solve  the  eigenvalue 
problems  for  any  parallel  or  nearly  parallel  axisymmetrical  flow.  Because  of 
some  simplifications,  we  may  divide  the  boundary  value  problems  into  two  dif¬ 
ferent  classes  as  follows: 

(i)  Flow  within  a  finite  domain.  In  this  case  both  the  end  points  |  ^ 
and  ^  are  finite.  One  of  the  two  limits  may  be  the  axis  of 

A 

symmetry,  i.e. ,  r  =  0,  or  ^  =  1.  The  boundary  conditions  are 
those  given  in  equation  (2 . 20)  which  may  be  expressed  in  terms 
of  ^  as 


(5.1) 


df 


The  factor  Vl  -  4  associated  with  ^  dropped  out  in 

the  determinantal  equation.  Hence,  we  may  use  the  boundary 
conditions  in  the  form  of  equation  (5.1).  It  should  be  noticed 
that  on  the  axis  of  symmetry,  the  boundary  conditions  are  still 


f  (1)  =  V  (1)  =  0 


because  the  limiting  conditions  that  the  radial  velocity  component 
on  the  axis  of  symmetry  vanishes,  give  these  two  conditions. 

This  is  one  of  the  essential  differences  of  the  axisymmetrical 
problems  from  those  of  symmetrical  flow  problem  of  the  two- 
dimensional  case. 

For  the  latter  case,  only  one  of  the  conditions  f  (0)  or  f*  (0) 
will  vanish  depending  whether  the  disturbance  is  symmetrical  or 
antisymmetrical. 
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The  determlnantal  equation  corresponding  to  the  boundary  conditions 
(5. 1)  is 


'll 

'21 

'31 

'41 

'l2 

^22 

^32 

'42 

II 

'll 

'll 

'll 

'll 

I 

'I2 

f 

22 

'I2 

=  0 


(5.2) 


where 


'll  -  'i<«i'  ■  'll  =  'l<?i».»*<=- 


2 

Similar  to  the  two-dimensional  case,  if  we  neglect  the  terms  of  the  order  of 
R^)  ^  and 


exponent 


-  f  ^  0  a  R  (W  -c)  d  i 


equation  (5.  2)  becomes 


(a^,  c) 


f' ,  "  (a,.  C) 


42 

f 

42 


(5.3) 


where 


f 

f 

f 

f 

11 

12 

11 

21 

— 

f 

f 

f’ 

f 

21 

22 

12 

22 

f  ’ 

f 

f  ' 

f  • 

11 

12 

11 

21 

f  ' 

f„ 

''4  = 

f 

f  ' 

21 

22 

12 

22 

(5.4) 
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Equation  (5. 4)  is  exactly  the  same  as  that  of  the  corresponding 
problem  for  two-dimensional  flow  between  solid  walls  in  relative 
motions.  * 


(ii)  Flow  with  an  infinite  domain.  It  is  the  case  for  wake,  jet,  or 
boundary  layer  flow  on  the  outside  of  a  circular  cylinder  or  a 
body  of  revolution.  In  this  case,  one  of  the  end  points  ^  or  ^ 
will  be  minus  infinity,  i.e. ,  r  =  “,  Co  =  Hence,  one  of  the 

viscous  solutions  0  '  ^  is  infinitely  large.  Our  boundary 

^  (0) 
conditions  require  that  f  should  be  a  linear  combination  of  f  ^  , 

f^^^^  and  alone,  i.  e. , 


f  4.  (0)  f  (0)  .  wi  (0) 

f  =  c  f  '  '  +  c  f  '  +  c  ' 
11  22  44 


(5.5) 


The  boundary  conditions  at  the  finite  limit  point  C  ^  "  are 


still 


f  (C^)  -  Vi  f  =  0 


(5. 6a) 


At  C  o  ~  ~  ”  I  the  fundamental  solution  0 


(0) 


tends  to  be  zero  very 


rapidly.  The  two  inviscid  solutions  behave  like  e 


ar 


The 


condition  that  f  tends  to  be  zero  as  r  tends  to  be  infinite  excludes 

-I-  r  -  r 

the  integral  e  Hence  f  must  be  proportional  to  e  for 

r  >  r.  .  This  condition  may  be  expressed  in  terms  of  ^  as  follows 
b 


V  1  -  ^  •  f  -«^f  =  0,  c. 


(5.6b) 


where  ^2  “  take  r^^  as  the  boundary  layer 

thickness  or  the  width  of  a  wake  or  a  jet. 


♦  Equation  (5.3)  corresponds  to  equation  (6. 13)  of  Reference  2. ,  but  there  seems 
a  printing  error  of  equation  (6. 13)  of  Reference  2,  i.e. ,  f^  should  be  f^.  Our 
G's  are  f's  in  Reference  2. 
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Ihe  determinantal  equation  corresponding  to  the  conditions  (5. 6) 
is 


f 

f  , 

f 

11 

21 

41 

f‘ 

fl. 

f 

11 

21 

41 

V^. 

f  -  a  f  \/i  -£  f  -  a  f 

12  1  12  ^  ^2  22  1  22 

0 

(5.7) 


or 


(5.8) 


where  6  =  0  for  two-dimensional  case  and  6  =  1  for  axisymmetrical 
case. 
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SECTION  6  -  EQUIVALENT  PRINCIPLE 


Since  both  the  fundamental  solutions  and  the  determinantal  equations  for 
the  axisymmetrical  flow  in  the  variable  4  are  similar  to  those  of  the  two-dimen¬ 
sional  flow  in  the  variable  y,  we  may  find  some  qualitative  relations  between 
these  two  cases. 

For  the  two-dimensional  flow,  the  second  derivative  of  the  mean  velocity 
2  2 

profile  d  W/dy  plays  an  important  role  in  the  stability  problem.  We  expect 

2  2 

that  the  second  derivative  d  W/d  ^  of  the  axisymmetrical  flow  plays  a  similar 
2  2 

role  as  d  W/dy  in  the  corresponding  stability  problem.  For  instance,  the  case 

2  2  2  2 
of  d  w/dy  =  0  should  correspond  to  the  case  of  d  W/d  ^  =  0.  In  the  two-di- 

2  2 

mensional  case,  d  W/dy  =  0  represents  the  plane  Couette  flow  case  which  is 

stable  for  infinitesimal  disturbances  at  all  Reynolds  numbers.  In  the  axi- 
2  2 

symmetrical  case,  d  W/d  ^  -  0  is  the  Poiseuille  flow  in  a  circular  pipe  which 

is  also  stable  for  infinitesimal  disturbances  at  all  Reynolds  numbers. 

For  the  two-dimensional  flow  problem,  if  the  basic  flow  has  a  point  of 
inflection  in  the  flow  field,  i.e. ,  W’*  (y  )  =  0  for  y  y  y^,  the  flow  is  very 
unstable. 

The  minimum  critical  Reynolds  number  of  the  basic  flow  with  a  point  of 
inflection  is  much  smaller  than  those  without  a  point  of  inflection.  For  instance, 
the  two-dimensional  wake  or  jet  flows  are  much  more  unstable  than  the  boundary 

layer  flow.  The  corresponding  case  in  axisymmetrical  flow  should  be  W’*  (4  )  =  0 

2  2  ^ 
for  ^  Even  though  d  W/dr  =  0  for  r  =  r  where  r  <  r  <  for 

^l^c  2  c  lc2 

an  axisymmetrical  wake  or  jet  d  W/d  4  for  ^  ^  .  For  instance,  for  the 

1  2 

far  downstream  wake  behind  a  body  of  revolution,  the  basic  velocity  profile  is 

W  (r)  =  1  -  b  exp  (-r^)  (6.1a) 

or 

W(C)  =  1  -  bexp  [-r/  (1  -o]  (6.1b) 
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where  b  is  a  constant  in  the  stability  analysis  (see  Part  II).  The  flow  field  lies 

between  r  *  0  and  r  =  “  or  ^  =1  and  ~  ~  " 

1^1  2 


dV 

dr^ 


=  2  b  (1  -  2r^)  exp  (-r^) 


(6.2) 


Hence, 


d^W 

dr^ 


=  0  when  r,  =  0  <  r  =  7=-  <  r_  =  ®  (6.3) 

1  c  V2  2 


2  4  r  2  T 

However,  d  W  =-br  exp  |_-r  (1-4)1  ^0 

9  C  C  ”* 

d4 


(6.4) 


when 


1  >  ^  >  -  00 


We  thus  expect  that  the  axisymmetrical  wake  or  jet  will  be  more  stable 
than  the  corresponding  two-dimensional  wake  or  jet  because  there  is  no  ’^effective” 
point  of  inflection  in  the  mean  velocity  profile  for  an  axisymmetrical  wake  or 
jet,  i.e. ,  d^W/d  4^  0  for  4  >  £  >  £2' 

The  axis  of  symmetry  in  the  axisymmetrical  flow  behaves  quite  different 
from  the  axis  of  symmetry  in  the  two-dimensional  flow  with  symmetrical  profile 
because  the  boundary  conditions  for  these  two  cases  are  different.  In  the  axi¬ 
symmetrical  case,  the  boundary  conditions  on  the  axis  of  symmetry  are  the  same 
as  those  on  a  solid  wall,  i.  e. ,  both  f  and  f*  are  zero.  For  the  two-dimensional 
case,  only  one  of  the  functions  f  and  f’  vanishes  on  the  axis  of  symmetry  depending 
on  whether  the  disturbance  is  symmetrical  or  antisymmetrical.  Hence,  the  axis 
of  symmetry  behaves  like  a  solid  wall  in  the  axisymmetrical  case  but  it  is  not  so 
in  the  two-dimensional  case.  As  a  result,  in  the  stability  analysis  of  an  axi¬ 
symmetrical  flow,  the  viscous  solution  always  takes  care  of  the  effect  of  vis¬ 
cosity  no  matter  whether  there  is  a  solid  wall  or  not.  We  may  draw  the  following 
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equivalent  principles  between  the  two-dimensional  and  the  axisymmetrical  cases 
of  stability  analysis: 

(a)  For  flow  of  finite  domain,  the  stability  problem  of  an  axisymmet¬ 
rical  flow  corresponds  to  the  stability  problem  of  a  two-dimen¬ 
sional  flow  of  unsym  metrical  profile  when  the  variable  is  used  for 
the  variable  y,  the  two-dimensional  variable.  The  stability  equa¬ 
tion  for  both  cases  is  equation  (5.3).  If  W  (  ^ )  =  W  (y)  we  would 
expect  that  the  critical  Reynolds  number  for  these  two  cases  are 
approximately  the  same. 

Because  of  the  above  principle,  it  is  expected  that  the  stability 
problem  of  the  Poiseuille  flow  in  a  circular  pipe  will  be  the  same 
as  that  of  plane  Couette  flow  because  W  ( ^ )  =  A  +  B  4  and 
W  (y)  =  A  +  By.  Actually  we  know  that  both  cases  are  stable  to 
infinitesimal  disturbance  for  all  Reynolds  numbers. 

It  should  be  pointed  out  that  Pretsch^  did  not  use  the  correct 
determinantal  equation  (5.3)  to  investigate  the  stability  of  a 
Poiseuille  flow.  He  arbitrarily  discarded  the  solution  f  which  is 
large  but  still  finite  on  the  axis  of  symmetry  ^  1.  He  used  then 
the  functions  f^,  f^  and  f^  in  the  stability  calculation.  He  had  to 
discard  one  of  boundary  conditions  on  the  axis  of  symmetry.  He 
used  f*  =  0  as  the  boundary  conditions  as  the  axis  of  symmetry. 
Finally  his  determinantal  equation  in  our  notations  is 


F  = 
3 


f 

f 

f 

11 

21 

41 

f 

f 

f 

11 

21 

41 

f' 

V 

f 

12 

22 

42 

4  4 


Go  =  0 
41  2 


or 


G„ 


4  (6.5) 

Pretsch  found  that  the  Poiseuille  flow  is  stable  according  to  equation  (6.5). 
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(b)  For  flow  of  infinite  domain,  the  stability  problem  of  an  axisym- 
metrical  flow  always  corresponds  to  the  stability  problem  of  a 
two-dimensional  boundary  layer  flow  when  the  variable  ^  is  used 
for  the  variable  y.  The  stability  equation  for  both  cases  is  equa¬ 
tion  (5.8).  If  W  ( ^  )  =  W  (y),  we  would  expect  that  the  critical 
Reynolds  numbers  for  these  cases  are  of  the  same  order  of  mag¬ 
nitude. 

This  result  is  especially  interesting  in  the  consideration  of  the 
stability  of  an  axisymmetrical  wake  or  jet.  We  expect  that  the 
minimum  critical  Reynolds  number  of  an  axisymmetrical  wake  or 
jet  would  be  of  the  same  order  of  magnitude  to  that  of  a  two-di¬ 
mensional  boundary  layer  flow  which  is  much  larger  than  that  for 
a  two-dimensional  wake  or  jet.  The  minimum  critical  Reynolds 
number  of  Blasius  boundary  layer  flow  based  on  the  displacement 

thickness  is  of  the  order  of  500.  For  boundary  layer  flows  with 
4 

pressure  gradient,  the  minimum  critical  Reynolds  number  varies 
from  100  to  1000.  On  the  other  hand,  the  minimum  critical 

5 

Reynolds  number  of  a  two-dimensional  wake  or  jet  based  on  the 
width  of  the  wake  or  jet  is  of  the  order  of  10.  Hence,  we  expect 
that  the  minimum  critical  Reynolds  number  of  an  axisymmetrical 
wake  or  jet  should  be  at  least  an  order  of  magnitude  higher  than 
that  of  a  two-dimensional  wake  or  jet. 

Numerical  calculations  will  be  made  to  check  the  above  principles.  It 

should  be  pointed  out  that  no  calculation  of  the  stability  of  a  two-dimensional 

boundary  layer  flow  of  velocity  profile  similar  to  equation  (6. 1)  has  been  carried 

out.  The  velocity  profile  (6, 1)  may  be  considered  as  a  uniform  velocity  (1-b) 

superimposed  on  a  variable  velocity  W.  ,  i.  e. , 

b 


or 


W  =  1  -b  +  W, 


=  b  W 


(6.6a) 
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Now  if  we  replace  c  by  c  such  that 


I 


be  =  c  -  (1  -b) 


(6.7) 


The  stability  equation  (4.  2)  holds  true  for  W  and  c  except  that  the  Reynolds 
number  should  be  replaced  by  R^^^  which  is 


R 


lb 


R^b 


(6.8) 


In  other  words,  if  we  use  the  actual  velocity  defect  bU  which  is  the  difference 
between  the  velocity  at  the  edge  of  the  wake  and  the  velocity  on  the  axis  of  sym¬ 
metry  as  our  characteristic  velocity,  we  may  use  the  velocity  profile  W  to  in¬ 
vestigate  the  stability  of  the  wake.  We  shall  use  this  transformation  in  our  in¬ 
vestigation  of  a  wake  flow  in  Part  II.  W  is  closer  to  the  ordinary  boundary  layer 
profile  than  W, 
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The  stability  eqiiation  of  an  axisymmetrical  flow  of 
an  incompressible  fluid  subjected  to  rotationally-symmet*^ 
ric  disturbances  is  derived  and  studied.  It  is  shown  that 
the  axisymmetrical  flow  is  always  stable  with  respect  to 
disturbances  of  the  tangential  velocity  component.  Hence 
jjthe  stability  equation  for  axisymmetrical  disturbance  is  ^ 
then  investigated  in  detail  for  arbitrary  basic  velocity 
profiles.  It  was  found  that  the  stability  of  an  axisymmet* 
rical  flow  may  be  reduced  approximately  to  the  corres¬ 
ponding  problem  of  two-dimensional  flow  by  proper 
transformations*  The  comparison  of  stability  of  an 
axisymmetrical  flow  with  that  of  a  corresponding  two- 
dimensional  flow  is  discussed  according  to  an  equivalent 
principle.  It  is  found  that  the  minimum  critical  Reynolds 
number  of  an  axisymmetrical  wake>  or  jet»  is  of  the  samel 
order  of  magnitude  as  that  of  a  two-dimensional  boundary 
layer  flow  which  is  much  higher  than  the  value  for  a  two- 
dimensional  wake,  or  jet. 
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